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Radiation pressure forces, the anomalous 

ACCELERATION, AND CENTER OF MASS MOTION FOR THE 
TOPEX/POSEIDON SPACECRAFT 

Daniel G. KubitschekW George H. Born* 

‘ Shortly after launch of the TOPEX/POSEIDON (T/P) spacecraft (s/ck foe 
Precision Orbit Determination (POD) Team at NASA’s Goddard Space Flight 
Center and the Center for Space Research at the University of Texas, discovered 
res idual alongtrack accelerations which were unexpected. Here, we describe the 
anal ysis of radiation pressure forces acting on foe T/P s/c for the purpose of 
; understanding and providing an explanation for foe anomalous accelerations. 

The radiation forces acting on foe T/P solar array, which experiences waxping 
| due to temperature gradients between foe front and back nur&ce^are analysed 
and die resul ting alongtrack accelerations are determined* Characteristics 
similar to those of the anomalous acceleration are seen* This analysis led lo the 
l development of a new radiation force model, which includes solar array warp ing 
and a solar array deployment deflection of as large as 2 °. As a result of this new 
model estimates of the empirical alongtrack acceleration are reduced in 
magnitude when compared to the GSFC tuned macromodel and are less 
* dependent upon the location of the sun relative to die orbit plane* If these 
r results are believed to reflect die actual orientation of the T/P solar array, then 
motion of the solar array must influence the location of the s/c center of muss. 
Preliminary estimates indicate that the center of mass can vary by as much as 3 
1 centimeters (cm) in the radial component of the s/c’s position due to rotation of 
the deflected, warped solar array panel The altimeter measurements rely upon 
accurate knowledge of the center of mass location relative to the s/c frame of 
reference. Any radial motion of the center of mass directly affects the altimeter 
l measurements. 


INTRODUCTION 

A few months after launch of the TOPEX/POSEIDON (T/P) spacecraft (s/c), the 
Precision Orbit Determination (POD) Team at NASA’s Goddard Space Flight Crater 
(GSFC) and the Center for Space Research at the Univeristy of Texas, observed residual 
alongtrack accelerations that were unexpected and of unknown origin . Figure 1 shows 
the residual alongtrack accelerations for both the GSFC prelaunch macromodel and me 
GSFC tiined macromodel (nonconservative force model was “tuned” using on-oibit 
observations up to and including cycle 16). These accelerations represent the residual 
constant ‘alongtrack accelerations, which are estimated in the orbit determination process 
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* . one. per day frequency. In addition, 

vary sinusoidally at a frequency of 1 cycle p er ThcscTOidual 

revolution or 1 CPR), were estimated m the o : J^[ v av{:ra g c constant 

accelerations have dm larges. Smtion” by 
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■r^S^eToS ' tpemnpt^ *^2 

^^o^'^frrtby adep^c upon bod. dm attitude Peering 
mode and/or the sun-orbit geometry. 


residual alongtrack acceleration 



Figure 1 Daily average 


f ft nT lant alongtrack accelerations for the first 160 days of the 
TOPEX/POSEIDON mission. 


T/P' orbits in one of four distinct yaw modes, depending upon the value of p. Table 
* i ’it. ... xu T x/P that is reauired to sunult&ncously m&intHin 

S‘£Ss5rs£*sss^s 

vector^ repr esents the angle between (he snn-line direction vector and the orbital plane, 
OkMM as L orbit angle teladve to the sunrise dirccnon veernr, and die 
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angle ¥ specifies the s/c yaw angle (rotation about the x* axis of the s/c frame of 
reference). 

Table 1 

TOPEX/POSEIDON YAW STEERING ALGORITHM 


Yaw Mode 

Orbit Forward Sinusoidal 
Orbit Forward Fixed 
Orbit Backward Fixed 
Orbit Backward Sinusoidal 


P' Region 
p'>+15° 
+15° £ p' > 0° 
0°>p'S-15° 
-15° >3' 


Yaw Angle 

- 90°+ (90° -- pOcos© 
T' = 0° 

180° 

= _90°- (90 c + PQcosQ 


In precise orbit determination of altimetric s/c, the radial component (or height 
compon^h) of the s/c’s location is the most important smce orbit 1 ™ 

this component directly affect the altimeter measurement. 

of the mission, two major efforts were launched to improve the predicted 1 ’ 

(1) development of an improved gravity field model 3 and (2) development of a detailed 
radiation force model for the T/P s/c 4 . Antreasian and ^Rosborougb (1992) 
analyzed a prelaunch micromodel approximation of the T/P s/c, which 
complex s/c geometry and attitude steering profile to determine thermal history 
SSion fof eachSirface. The acceleration profiles due to incident *£ patron 
pressure, (Earth albedo and infrared emission, and thermal emission of radiation from 
lach surface of the s/c were computed for various values of P'This munmodd was too 
computationally intensive for operational orbit determination, but provided the necessary 

information needed for development of a much simplified radiation force model. 

The POD Team used the accelerations from the micromodel analysis to den 
macromodel approximation known as the “box-wing” model. This of 

eight surfaces, each having a set of property parameters which were adjiisled so that the 
resulting accelerations matched the micromodel accelerations . Figure 2 shows th 

evolution of the T/P nonconservative force model. , 

Upoii observation of the anomalous acceleration, deficiencies in the macromodel 
were immediately suspected. Solar array warping was known to exist prior tolaunclL bu 
thought to be negUgible. In order to gain some insight it is important to und«stand the 
* orientetili of thfs/cduring each of the yaw steering modes. During the orbitforward 
fixed yaw mode the s/c remains fixed relative to the orbit frame of reference (defined to 
be the local vertical local horizontal frame) with the +x*fc axis pointing m the 
the local horizontal direction vector (the projection of the velocity vector onto th 
horizontal direction vector is positive). This orientation represents a yaw ang e o • 
During the orbit backward fixed yaw mode, the yaw angle remains fixed at 180^ The 
siSdil yaw modes are more complicated. During toe orbit forwato 
mode, the -y s/c axis of the spacecraft oscillates back and forth *n:oss the direction 
local horizontal while the amplitude of oscillation vanes as a function of P . The orbit 
y»w ml has the osdUmmg *»«. 

direction vector. Here, warping of the solar array or a deploymentdefletf^ 
array could be responsible for residual alongtrack accelerations. If it is assumed that th 
solar array is flatten only the area associated with its thickness would lie projected in 
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acceleration of 0.2 - 0.4 nm/a’ (1 om - 1 nnom«er - 1*10* ™w») 

to acoontmodate the anonatlon accdentjon, the POD T«m hhj«i 
(adjusted) the surface property parameters to decrease the root-mean - siiueu ) 

difference between die computed lie locations and the observed (via tractate &») 
location This was expected to remove most of the mismodelmg associated 
simplified macromodel approximation. In addition, constant alongtrack »o^crati^ ? 
HTmn/s 2 and 0 20 nm/s 2 in the and +** directions were estimated from the 
residual alongtrack accelerations and applied throughout the entire 160 day 
toe timed macromodel results shown in figure 1. This does mdeed i^veafew oi 
discontinuities seen at yaw mode transition, but does not change the characteristics seen 

oniy group » ob^o an .«*. , bdmvjor. TheT/P 
Navigation (NAV) Team at NASA’s Jet Propulsion Laboratory (JPL) also obswved 
ZS Semifixed accelerations 7 . TTte NAV Team uses a 

i Tc^ noXcount for variations in surface properties « motion of toe actor 
the anomalous satellite-fixed accelerations discussed by Frauenholze/ al. (1993) are ery 
different from and much larger than the accelerations discussed in this paper . 

- t 

RADIATION FORCE CALCULATION 

The 'computational model for solar array warping consists of a flat p!ate, witti the 
properties shown in table 2, which is decomposed into 100 partitions. Nodes > are Jucat 
at toe center of each partition. Warping of toe solar array results from temperatoe 
differences between toe front and back surfaces of the solar array and this war P^ g ^V , 
quantified using the physical dimensions and the effective ^efficient of thermal 
expansion for toe composite solar array*. Curvature along toe lengto ' 
width of the panel may then be determined for any given temperature difference via the 

following equation: 


R = 


a cte AT 


Figure 3 shows the compound curvature of the solar array for a t^perature dffoen^ of 
10 K. The resulting tip deflection is approximately 20 cm. For toe analysis presented m 
this pamir only curvature along toe lengto of toe solar array is cohered 

Several MATLAB routines were developed to simulate toe T/P orbit, s/c onentanon, 
and compute accelerations due to radiation pressure usmg s/c position, floaty, 
temperature, and attitude telemetry which are extracted from toe Selected Telemetry 
Record^ (STR) files available at JPL. These data are also used to define toe unportant 
referencte frames, determine Earth shadow exit and entry, 
radiation forces acting on toe warped solar array are computed ; ^ 
accelerations are determined in the orbit frame of reference. The ^difference jbrtwmi toe 
warped Uar array accelerations and toe flat solar array accelerations is detwmmcd and 
* - • . i ixi- 4 ., nflMHtmfinn thflt thfi solar array is flat, ine 


warped Solar array accelerations anu me um ^ — ; ; fl 

represents toe error associated with the assumption that toe solar array is flat. 
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difference in the alongtrack component was averaged daily over a period covering a p 

lay* of *• operational «HL these leys were chos« based 

upon STR availability at the time. 


Table 2 


TOPEX/POSEIDON SOLAR ARRAY PHYSICAL CHARACTERIS1 ICS 


Vihie 


Definition 

Length 

L 

7.72 m 
3.30m 

Width 

w 

A 

0.0349 m 

Thickness 

a 

U 

2300 kg 

s/c Mass 


23.4x10'* K' 1 

Coefficient of Thermal 

&cte 


Expansion 

Ct 

q.go 

Absorptivity 

V 

0,20 

Reflectivity 

r 

0.20 

Specularity 

) 

p 

0.81 

Cell SideEmissivity 

e / 


Back Side Emissivity 


(1,85 



i 


Figure 3 Computational model for solar array warping simulations. 


Figure 4 shows the daily average constant alongtraok acceleration due to sole ary 
warping. Examination of the temperature telemetry shows a difference of V 

10 K - 12 K between the front and back surfaces of the solar array during p 
munition: According to the watping model used in this m+m. a plane winch ts 
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oriented tangent to the center point of the aote «ay. fJT^ge 

** -»** “ — — •“ *• 

results obtained by Kar and Ries (1993). 

daily a vbraob alongtrack acceleration 



Figure" -t 


Predicted daily average constant alongtrack accelerations doe to solar aiTay warping. 


Comparison of the daily average alongtrack acceleration due to solar 
with the /anomalous acceleration results displayed in figure 1 show s ® r , c , 

during sinusoidal yaw modes. The most promising results are shown during 5 *e second 
orbit forward sinusoidal yaw period (days 290 - 340) where the daily averag gtr 
accelerations remain relatively constant During that period p' caches a maximum ue 
of approximately 50° (never experiencing fall sun orbits). TheP values ahownmfipg^ 
1, Xch corresponds to the orbit forward sinusoidal yaw period, ate do not exceed 55.7 
and the associated residual alongtrack accelerations exhibit a constant trend. 

aS atmparison of die results obtained from this solar array 
shows rimilar characteristics throughout the sinusoidal yaw mode regions it tannot 
observed during the fixed yaw periods. Also the mapnfadeof fae 
Sgtrock accelerations due to solar array warping area taxor of Am 
of&ialous acceleration. Nonetheless, these resuite provide j^fication enough 
proceed' with modification of the macromodel being used by die POD Team. 

MODIFICATION of the topex/poseidon macromodel 
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The original priori and toned 

reconstruct foe published immedtotoly notod 

relative to which changes may be evaluated. p • vn i aunc Vi the 

and corrected. The first had to do with the prelaunch mac^odeL pn ™* m 2 
T/P solar array design specifications listed the solar array 2 

w^ changed just before launch to the ^ surfa^^ of^ 
The second discreoancy had to do with the tuned macromodel. The tuning prooas is 
empirical in nature and can result in unrealistic changes to 

values for two of the surface property parameters adjusted dunng g pr 

, Glided that S four panel awroximation ° f 

1 used in tii modified macromodel approximation as seen tn figure 5. A few gumptions 
Jave been made in order to simplify the solar array warping model. (1) the 
composed of 5052 aluminium honeycomb material, (2) the value for the 
Slpansion is the value reported in table 2 (3) the 

the solar arrav is neglected, and (4) each of the solar array jomts are folly extended and 
residual alongtrack accelerations were determined over the 

oneratioiial orbit Figure 6 indicates that solar array warping alone cannot folly explain 
th^omSous acceleration during the sinusoidal yaw modes, torn 
consistent with foe results obtained in sin.ulat.ons .re » "to SHS 
during foe fixed yaw periods are due to adjustments made to foe X+ and X- emissivity 

«ny deployment angle emits are another possibility stao. d^oyront of foe 
solar array could not be guaranteed to better than ±1 . After tabbing 
engineers it was agreed that a 2" deployment deflection angle would not be unrealistic^ 
n macromodel was further modified to include a 2" defleet.cn angle m 
SfoSo solar .harping, to other words, foe entire structure is routed about a 

lateral axis (axis orietoedlmross foe width of foe soler array) attached to the aid of foe 
S X wtocVlXeted nearest to foe s/ebus. This deflecrion is assumed to ronam 

constanfat all times. The resulting constant 

in fioure 7 and indicate that most of foe remaining CAAs dunng foe sinusoidal ya 
SSiJSJSTfcr. Howeva, a negative meat CAA ranains throughout all yaw 
modes for which no satisfactory explanation has yet been found. 

EVALUATION OF THE MODIFIED MACROMODEL 
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Figure 5 The four 

i 


panel approximation of solar array warping under a 12 K tempera ture difference 
between the front and back •orfacei. 



i 

Figure 6 The CAA for solar array warping, X+ emisiivity « 0.3, and X- emuaivity - 0.9. 
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CONSTANT ALONGTRACK ACCELERATION 



too wo 

Days In Operational Ortrit 

Figure 7 The CAA for solar array warping, a solar array deflection angle of 2% X+ emissivtty - 0J, 

* ( and X- emissivity - 0.9. 

for both the SLR (satellite laser ranging) and DORIS (Doppler OrWtopaphy toge 
Integrated System) observations show close agreement for each 
™SLls» T is important to point otit that we wotdd not expect the nae of an 

improved macromodel to result in improwrf fits to the 

enSrical accelerations in all three cases. These empirical accderanon pai ameters ad to 

absorb tmmodeled accelerations. However, we can evaluate the ’° f 

j„i u v m otead lookin g at die magnitude of the empirical accelei ahons. 

rrtmSt^d, dm R^fita 

While the magnitude of the CAA parameter is greatly improved through a reducb 

ma8 C!^tSrS te C SSof the modified 

SSThave Sexumded through cycle 52. Tdtle 4 shot™ 

! 1-52 forthe GSFC tuned macromodel, the modified maciomodel which has a solar jwray 
deflection angle of 2°, and the modified macromodel which has solar array deflection 

i of;l.5°and SA+ 

l^ C is“s" £ £ ££• macromodd whid, ossa a so.ar 
array deflection of 1 *5°. 
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The dependence of the CAA values on p' were also examined. Figure 8 shows the 
CAA values as a function of p' and demonstrates that the modified macromodel removes 

much of the p / dependence. 



P', degrees 


Figure 8 The CAA as a function of P' for three different macromodels during cycle 1-52. 


11 




01/04/02 11:25 FAX 818 393 6388 


JPL NAV SYS 


S] 014 


Table 4 

CONSTANT ALONGTRACK ACCELERATION STATISTICS FOR THREE DIFFERENT 
MACROMODELS OVER CYCLES 1-52. 


Macromodel 

Mean 

(mn/s 2 ) 

Min 
( mn/s 2 ) 

Max 

(cm/* 2 ) 

Standiird Deviation 

(nm/» z ) _ 

GSFC Tuned 

-0.158 

-0.875 

0.549 

0.217 

Modified 2.0° 

-0.042 

-0.557 

0.466 

0.195 

Modified 1.5° 

-0.036 

-0.625 

0.321 

0.159 


MOTION OF T HE TOPEX/POSEEDON CENTER OF MASS 

The results shown thus far laid support to the hypothesis that warping and a 
deployment deflection of the T/P solar array may be responsible for the anomalous 
acceleration. If this is indeed the case, then we must take into consideration motion ot 
the s/c center of mass. The T/P solar array contributes approximately 10% to the entire 
mass of the s/c. Solar array deflection and solar array warping results m a tip deflection 
of almost 24 cm. Deflections at die location of the solar array center of mass can be as 
much as 15 cm which would displace the total s/c center of mass by approximately 1.5 
an. This displacement may be determined as a function of time provided that knowledge 
of the T/P solar array pitch angle and solar array front and back surface temperatures are 
available. The MATLAB routines used in the constant alongtrack acceleration 
simulations were modified to compute die s/c center of mass displacement (or offset) for 
three characteristic p' regions: (1) high p' during which the s/c is illuminated throughout 
it’s entire orbit, (2) moderate p' during which the s/c is still in the sinusoidal yaw mode, 
but now it experiences periods of Earth occultation during each orbit, and (3) low P' 
where the s/c is in the fixed yaw mode with the solar array panel rotating at a frequency 
of 1 per satellite revolution. Figure 9 shows the displacement of the T/P colter of mass in 
the radial direction. The radial direction is the most important since displacements in the 
radial component directly affect the altimeter measurements; it is assumed that the cento 
of mass! location is known relative to the altimeter measurement reference datum. These 
computations show that the center of mass can vary by as much as 3 cm in the radial 
component during fixed yaw. During sinusoidal yaw, the displacement tends to be more 
characteristic of a negative constant bias. The center of mass of&et is determined relative 
to the s/c frame of reference, so a negative displacement moves the coiter of mass away 
from t^e surface of the Earth. Neglecting this displacement results in altimeter 
measurements that are short, thus placing the sea surface higher relative to the center of 

file Earth. ... .... . 

The difficulty now lies in observing the coiter of mass motion if it exists. Altimeter 

bias estimates at NASA’s primary verification site (Harvest platform) woe corrected 
point by point for the suspected center of mass displacement thinking that this might 
reduce the standard deviation of the bias estimates. No reduction was observed, however, 
a shift in file bias estimates of approximately 1-2 an was seen. This leads to an inclusive 
result, but provides another example of the importance of understandin g this motion and 
how it affects the altimeter measurements if at all. Using the GEODYN II attitude and 
temperature model, a time series of the center of mass offset has been computed for 
cycles 5-121. Figure 10 shows the average, by cycle, center of mass offset for cycles l- 
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100. The displacement tends to be negative and correlated with P' as expected. 
GEODYN II was then modified to accept these offsets for cycle 46 and apply the onset 
to correct for the location of the center of mass. This led to a degradation in thefitsto the 
tracking data. The SLR RMS of fit changed from 0.0275 m to 0.0283 m. The DORIS 
RMS of fit was unchanged at 0.0562 cm/s. The magnitudes of the SLR RMS of fits are 
approximately 2 cm larger than the current state of T/P fits to the tracking data. This is 
due to the use of early GEODYN setups known as the REPRO-1 setups. 

Using the time series information, the center of mass offsets were tied to the T/P 
reference groundtrack latitude and longitude in an effort to investigate the possibility of a 
geographic correlation. This data was then transformed to a gridded data set and 
interpolated to produce a surface representing the error induced by motion of the center 
of mass. A global map for the each of the first 20 cycles was created and shows features 
that are latitude dependent during periods of high P' as seen in figures 1 1 and 12. During 
the orbit fixed yaw modes there are no latitude dependent features and the center of mass 
offset appears to be nothing more than a source of noise as seen in figure 13. Figure 14 
displays the value of p' over the first 20 cycles. 


i 



Figure 9 The Radial motion of the TOPEX/POSEIDON center of mats over one orbit revolution due 
to a 2° solar array deployment deflection and due to tolar array warping. 
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CM offset, oa 


Figure 10 Average cycle by cycle center of mass offset for cycles 1-100. 



Figure 11 Geographic surface plot for cycle 1 derived from the center of mass; offset time series. 




CM offset, an 
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CYCLE 017 



Figure 12 Geographic surface plot for cycle 17. 


CYCLE 007 



Figure 13 Geographic surface plot for cycle 7. 
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Figure 14 f}' for cycles 1-20. 


CONCLUSIONS 

t 

The original purpose of the research discussed in this paper was to investigate the 
anomalous acceleration and provide and explanation for these accelerations. The theory 
that solar array warping was responsible for die residual constant alongtrack accelerations 
during fhe sinusoidal yaw modes was investigated and the results shower! that solar array 
warping alone cannot account for all the residual accelerations. Solar array deployment 
deflections were then studied showing that deflections of approximately 1.5° - 2° could 
account for die remaining residual constant alongtrack accelerations during die sinusoidal 
yaw modes. The anomalous acceleration seen during die fixed yaw modes was removed 
by adjusting die X+ and X- surface emissivities in the macromodel. It is important to 
look at the performance of the modified macromodel during later cycles to be sure that 
the resulting accelerations are improved over a large number of observations. The plot of 
the estimated constant alongtrack accelerations, shown in figure 8, demonstrate the 
ability of the modified macromodel to remove much of the p' dependence. However, if 
these Results are to be believed, then we must look at the effects of solar array motion 
upon the s/c center pf mass. Several attempts have been made to observe this motion, 
m tracking data residuals degrade slighdy when the center of mass offsets are included, 
(2) bids estimates at Harvest platform do not show a clear reduction m standard deviation 
of the estimates over several cycles, hut the center of mass motion would be observed as 
a systematic error in the altimeter calibrations. Cycle by cycle surface plots, which 
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represent the center of mass offsets geographically, show some interesting features. 
These features are on the order of 1-2 cm and may not be observable in the altimetry. It 
may turn out that this motion has no significant affect upon the T/P altimetry system, but 
until all possibilities have been exhausted, the question cannot be laid to rest. 
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